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' Abstract. In a consistently formulated pQCD framework incorporating non- 

zero mass heavy quark partons, there is still the freedom to define parton distri- 
^ ' butions obeying either mass- independent or mass-dependent evolution equations, 

0> , contrary to statements made in a recent paper by MRRS. With properly matched 

'/^ ' hard cross-sections, different choices merely correspond to different factorization 

schemes, and they yield the same physical cross-sections. We demonstrate this 
principle in a concrete order calculation of the DIS charm structure function. 
1^ , We also examine the proper matching between parton definitions and subtrac- 

0^ ' tions in the hard cross-section near threshold where the calculation is particu- 

larly sensitive to mass effects of the heavy quark. The results obtained from the 
Q_i' general-mass formalism are quite stable against different choices of scale and ex- 

hibit a smooth transition in the threshold region (using either mass-independent 
or mass-dependent evolution), in contrast to results of another recently proposed 



r~| , scheme. 

> : 

^ ■ Recent improved measurements of heavy quark production in both leptopro- 

■ duction and hadroproduction demand better theoretical understanding of the 

underlying QCD physics. One must consider the changing role of the heavy 
quark (denoted by H) from threshold to the very high energy limit available 
in current and future experiments. Specifically, we address the question of 
the mass dependence in the evolution of the parton distribution functions 
(PDF's), /„/p(x,/x). 

In the ACOT factorization scheme, [1,2] the heavy quark mass (niH) is 
fully incorporated in the hard cross-section and threshold (/i = mn) matching 
conditions on and fa/p{x,ii). This scheme has the desirable feature 

that the evolution of all partons, including the heavy quark, is controlled by 
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the mass-independent MS kernels; i.e., fa/p{x,ii) satisfy the well established 
evolution equations. This is a precise consequence of the renormalization 
scheme in the theory with mass, it is not an approximation.^ 

A general mass formalism can, in principle, be implemented with factoriza- 
tion prescriptions different from that of ACOT. In a consistently formulated 
theory, physical predictions from different schemes must remain the same up 
to higher-order corrections. In a recent paper, MRRS [3] proposed a procedure 
within the general mass approach involving mass- dependent splitting functions 
for the evolution equations. They questioned the validity of the ACOT scheme 
because of its use of mass-independent evolution. It is important to clarify 
this question and determine what schemes are truly viable, since further de- 
velopment of heavy quark theory and phenomenology rely on the availability 
of dependable calculations. 



I SCHEME INDEPENDENCE OF PHYSICAL 

CROSS SECTION 

We consider charm production in DIS in a general renormaliza- 
tion/factorization scheme R. The dominant contributions are:^ 

= ""fc/p ® ^SLc + ""fg/p ® ( 4r-cc- - ""a^J^Lcc-) + o{al) (1) 

^a"^ denotes the hard cross section in a given renormalization scheme R at 
order n. cr" denotes the scheme independent partonic cross section.'^ ^cr" de- 
notes the subtraction terms which removes the potential heavy quark coUinear 

singularities from cr". The physical cross section (Jphys is scheme independent. 
The choice of subtraction ^a" defines the scheme, and hence the parton distri- 
butions ^fa/p and their associated splitting functions ^Pa^b- In fact, to this 
order these terms are connected by the relation: ^cr^^^cc = ^/i/] ® ^l^*-^^ 
where ^/c/] is the 0{a\) perturbative distribution of c in g. For the complete 
(7 SUB we nave: 

^) As the MS scheme yields mass-independent renormalization constants to all orders, the 
evolution kernels must be mass-independent. As mH — > 0, the massive hard scattering ct's 
reduce to precisely the zero-mass MS cr's without any finite renormalization; in this sense, 
the ACOT scheme is uniquely the minimal massive extension of the MS scheme. [1] 

We refer to the term in Eq. 1 as: (Jtot = (^he + {'^hc — (^SUb), where axoT represents 
the total physical cross section, anE the heavy- flavor excitation term (7*0 — »■ c), anc the 
heavy-flavor creation term (7*51 — > cc), and asuB the subtraction term. 

To focus on the heavy quark parton singularities (of the form \n{mH/Q)), throughout 
the paper it is understood that the collinear singularities associated with light partons are 
always subtracted from cr" according to the usual MS scheme. 
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FIGURE 1. F2 vs. n for DIS c-production. a) ^3^^, F^^^ and the difference 
F^^ — F2^^ . The sohd curves are for the mass-independent evolution scheme, and the 
dashed curves are for the mass-dependent evolution scheme, b) F2^^ and F^^ — F^f^^. 
The difference between the mass-independent evolution and mass-dependent evolution for 
pTOT j^igj^gp order and comparable or less than the /^-variation. 

asuB = %/P ® ^aj;L., = %/p ® ^ rl ^ ^P«. ® (2) 
Near threshold (/i ~ m^y), we have 

auE = ""Up ® «aSLe - ^/./p ® ^ f ^ ^^iH. ® al°L, + O^al) (3) 

which reflects the fact that osub — cthe, by construction. This yields 
aroT — o'hc + C'(ag) independent of the specific choice of scheme R, and 
therefore independent of the choice of ^al^^^cc ^md ^P^^^. This is the key 
mechanism that compensates the different effects of the mass-independent vs. 
mass-dependent evolution, and yields a cttot which is identical up to higher- 
order terms. 



II MASS-DEPENDENT OR MASS-INDEPENDENT 

EVOLUTION? 

Ref. [3] strongly advocated the use of mass-dependent evolution of the par- 
ton distributions, and questioned the correctness of using the familiar MS 
mass-independent evolution in ref. [1]. Is this criticism justified? No-provided 
the proper matching between the hard scattering cross sections and the parton 
distributions is observed. Although the ACOT scheme chose to use a mass- 
independent evolution for its simplicity and natural connection to the familiar 
zero-mass parton results, this choice is not an approximation — the physical 
predictions are insensitive to this choice. We explicitly demonstrate this by 
implementing a mass-dependent evolution scheme in the same formalism, and 
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FIGURE 2. F2 vs. Q2 f^j. j^jg c-production. The scale choice is: n'^ = Q'^/A. a) 
The general mass formalism with mass-dependent evolution, b) The MRRS scheme with 
mass-dependent evolution. This figure is taken directly from ref. [3], Fig. 6b. The curves 
are re-labeled to correspond to our notation. 

compare the new results with the existing mass-independent ACOT calcula- 
tion. In the next section, we shall compare this new massive calculation with 
that of MRRS. 

In Fig. 1 we display the separate contributions to F2 for both mass- 
independent and mass-dependent evolution. The matching properties dis- 
cussed earlier are best examined by comparing the (scheme-dependent) F.j^^ 
and F2^^ contributions of Fig. la. We observe the following, i) Within each 
scheme, and F2^^ are well matched near threshold, ii) The matching of 
F^^ and F2^^ ensures that the scheme dependence of F^*^-^ is properly of 
higher-order in a^. 

The lesson is clear: the choice of a mass- independent MS or a mass- 
dependent (non-MS) evolution is purely a choice of scheme, and becomes 
simply a matter of convenience-t/iere is no physically new information gained 
from the mass- dependent evolution. 

Ill ^2™^ SHOULD MATCH FIXED-ORDER 
RESULT NEAR THRESHOLD 

A necessary consistency check of any proposed general mass scheme in the 
important intermediate energy region is that it yield the correct limits for 
threshold and asymptotic energies. Although this feature is, in principle, 
built into the formalism discussed in Eq. 2, it is essential to verify that the ex- 
pected behavior actually appears in specific implementation of the scheme — in 
particular, near the threshold region where the required cancellation between 
the F2^ and F2^^ terms is delicate 

Fig. 2 compares F2 vs. the physical variable Q using the same charm mass, 
scale choice, and format of an existing MRRS plot, (c/. Fig.6b, Ref. [3]). 
In Fig. 2a, F^'^'^ smoothly interpolates between F.^'^ in the threshold region, 
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FIGURE 3. a) Same as Fig. 2, but with /i^ = Q^- b) /^-variation in the general mass 
formahsm with mass-dependent evolution. The variation of F2'~'^ is small in contrast with 
that of F2^- Note, we use a linear scale for F2 to give a proper perspective. 



and in the high energy region, demonstrating the desired matching.^ In 
contrast, the corresponding curve in Fig. 2b (taken from Ref. [3]) shows a sud- 
den rise at threshold which indicates a lack of expected cancellation between 
F^^ and ¥2^^ just above threshold. Since both calculations use the same 
(MRRS) evolution kernel, this difference can only be due to the implemen- 
tation of the subtraction procedure: Fig. 2a uses the general mass formalism 
described above, [1,2]; Fig.2b uses the procedure described in Ref. [3]. 



IV Fi^^ SHOULD BE STABLE UNDER 
SCHEME/SCALE CHOICE 

Irrespective of smoothness of the theory calculation, a physical prediction 
must be stable against the choice of artificial theoretical parameters such as the 
renormalization and factorization scale Fig.3a shows our calculation (again, 
using the mass-dependent evolution) with a different choice of compared to 
Fig. 2a. In Fig.3b we note that while changed dramatically with choice 
of /i, this change is compensated by F2^^ such that F™^ remains stable. 

The corresponding comparison is not available for the MRRS procedure. If 
the threshold behavior evident in Fig. 2b persists for other choices of scale, the 
position of the sudden rise in cross-section will move with the choice of scale; 
that would be cause for concern.^ 



The intuitive reason for these physical limits is identical to the original mass-independent 
scheme. [1] In the threshold region the "heavy quark" behaves more like an extrinsic heavy 
object, while for high energies the "heavy quark" behaves like an intrinsic parton. 

For a scale choice of /i^ — Q^/k^, the charm quark evolution turns on at = k^rn^, 
where k^ is a arbitrary parameter. 
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V SUMMARY 



In conclusion, we find that the presence or absence of heavy quark mass 
in the evolution kernels has no physical efi^ect on cross sections — there is no 
physically new information gained from a mass-dependent evolution. 

A mass-dependent evolution scheme has inherent difficulties. The definition 
of the mass-dependent splitting kernels introduces a kinematic ambiguity since 
it is impossible to have an on-shell masslcss particle split into two collincar 
massive on-shell particles (e.g., g cc); this leads to a number of different 
definitions for the mass-dependent sphtting kernels. [5] Since the PDF's in the 
mass-dependent scheme are not MS, the hard scattering cross sections a must 
be converted from the MS scheme to this scheme. At NLO, this conversion, as 
well as the PDF evolution, requires two-loop splitting kernels (P*^^-') computed 
in the mass-dependent scheme; these results do not yet exist. 

The hallmark of any consistently formulated pQCD theory is: i) proper 
matching between F^^ and F2^^ at threshold, and ii) insensitivity of the 
physical F^"'^"^ to the arbitrary scheme/scale choice. The properly imple- 
mented general mass formalism satisfies these criteria. The ACOT scheme 
(with mass-independent MS evolution) has the added advantage that the 
known mass-independent MS kernels follow naturally from the renormaliza- 
tion scheme adopted. 
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